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Summary. — We discuss the influence of the cosmological constant A on the grav- 
itational equations of motion of bodies with arbitrary masses and eventually solve 
the two-body problem. Observational constraints are derived from measurements 
of the periastron advance in stellar systems, in particular binary pulsars and the 
solar system. For the latter we consider also the change in the mean motion due to 
A. Up to now, Earth and Mars data give the best constraint, A ~ 10~^^ km~^. 
If properly accounting for the gravito- magnetic effect, this upper limit on A could 
greatly improve in the near future thanks to new data from planned or already oper- 
ating space-missions. Dark matter or modifications of the Newtonian inverse-square 
law in the solar system are discussed as well. Variations in the behavior are 
considered in the form of either a possible Yukawa-like interaction or a modification 
of gravity of MOND type. 

PACS 04.80.Cc - 95.10.Ce, 95.30.Sf, 96.36.-hx. 



1. Introduction 

The understanding of the cosmological constant A is one of the most outstanding topic 
in theoretical physics. On the observational side, the cosmological constant is motivated 
only by large scale structure observations as a possible choice for the dark energy. In 
fact, when fixed to the very small value of ~ 10~^^ km~^. A, together with dark matter, 
can explain the whole bulk of evidence from cosmological investigations. In principle, 
the cosmological constant should take part in phenomena on every physical scale but 
due to its very small size, a local independent detection of its existence is still lacking. 
Measuring local effects of A would be a fundamental confirmation and would shed light 
on its still debated nature, so it is worthwhile to investigate A at any level. 

Up till now, no convincing method for constraining A in an Earth's laboratory has 
been proposed pj. Astronomical phenomena seem to be more promising. The cosmo- 
logical constant can affect celestial mechanics and some imprints of A can influence the 
motion of massive bodies. In particular, the effect on the perihelion precession of solar 
system planets has been considered to limit the cosmological constant to A ~ 10^"^^ km~^ 



© Societa Italiana di Fisica 



1 



2 



PHILIPPE JETZER AND MAURO SERENO 



So far local physical consequences of the existence of a cosmological constant were 
investigated studying the motion of test bodies in the gravitational field of a very large 
mass. This one-body problem can be properly considered in the framework of the spheri- 
cally symmetric Schwarzschild vacuum solution with a cosmological constant, also known 
as Schwarzschild-de Sitter or Kottler space-time. The rotation of the central source can 
also be accounted for using the so-called Kerr-de Sitter space-time. In [2] we carried out 
an analysis of the gravitational TV-body problem with arbitrary masses in the weak field 
limit with a cosmological constant. This study was motivated by the more and more 
central role of binary pulsars in testing gravitational and relativistic effects. 

Gravitational inverse-square law and its relativistic generalization have passed signif- 
icant tests on very different length- and time-scales. Precision tests from laboratory and 
from measurements in the solar system and binary pulsars provide a quite impressive 
body of evidence, considering the extrapolation from the empirical basis [31 3] . First 
incongruences seem to show up only on galactic scales with the observed discrepancy 
between the Newtonian dynamical mass and the directly observable luminous mass and 
they are still in order for even larger gravitational systems. Two obvious explanations 
have been proposed: either large quantities of unseen dark matter (DM) dominate the 
dynamics of large systems or gravity is not described by Newtonian theory on every scale. 

High precision solar system tests could provide model independent constraints on pos- 
sible modifications of Newtonian gravity. The solar system is the larger one with very 
well known mass distribution and can offer tight confirmations of Newtonian gravity 
and general relativity. Any deviation emerging from classical tests would give unique 
information either on dark matter and its supposed existence or on the nature of the de- 
viation from the inverse-square law. The orbital motion of solar-system planets has been 
determined with higher and higher accuracy 5J and recent data allow to put interesting 
limits on very subtle effects, such as that of a non null cosmological constant [21 [6l [7l [8] . 

2. The two-body problem 

The total Lagrangian for two particles can be written as 

(1) £ ~ ^TOa^a + G ^""^^ + ^rribvl + (5£pN(A=o) + SCa 

where x = Xa_ — Xh is the separation vector and (5>CpN(A=o) and SCa are the pN and 
A-contributions, respectively, with 

(2) 6Ca = ^{nii^xl + rribxl). 



Due to cosmological constant, the energy of the system is modified by a contribution 
—SCa- The pN and A corrections are additive and can be treated separately. Since the 
perturbation due to A is radial, the orbital angular momentum is conserved and the orbit 
is planar. The main effect of A on the orbital motion is a precession of the pericentre. 
One gets for the contribution to the precession angular velocity due to A, 

(3) a;A = ^^Vl-e2, 

where e is the eccentricity and Pb the Keplerian period of the unperturbed orbit. This 
contribution should be considered together with the post-Newtonian periastron advance, 



SOLAR SYSTEM TESTS OF THE COSMOLOGICAL CONSTANT 3 

Table I. - Limits on the cosmological constant due to extra-precession of the inner planets of 
the solar system. 

Name Sijj (arcsec/year) i^a (deg/ycar) Aum (km^^) 

Mercury -0.36(50) x lO"* 9.61x 10^^A/(1 km-^) 4x10-^'' 

Venus 0.53(30) X 10"^ 2.51xlO^'*A/(l km-^) 9x10"^^ 

Earth -0.2(4) x 10"^ 4.08x 10^'^A/(1 km-^) IxlO"^*^ 

Mars 0.1(5) X 10"^ 7.64x 10^''A/(1 km-^) 2x10"^^ 



cOpN 3(27r/Pb)^/^(GAf/c^)2/^(l - e^)"^ The ratio between these two contributions 
can be written as, 

(4) ^ _^Pa ^ 1— A 

where R — a(l — e^)'^/® is a typical orbital radius with a the semi-major radius of the 
unperturbed orbit, p = Mj {^-kB? j'S) is a typical density of the system and pA = c^K/^-kG 
is the energy density associated to the cosmological constant. The effect of A can be 
significant for very wide systems with a very small mass. 

2T. Interplanetary measures. - Precessions of the perihelia of the solar system planets 
have provided the most sensitive local tests for a cosmological constant so far [HI [TOj ■ 
Estimates of the anomalous perihelion advance were recently determined for Mercury, 
Earth and Mars [111 [S] . Such ephemeris were constructed integrating the equation of 
motion for all planets, the Sun, the Moon and largest asteroids and including rotations 
of the Earth and of the Moon, perturbations from the solar quadrupole mass moment 
and asteroid ring in the ecliptic plane. Extra-corrections to the known general relativistic 
predictions can be interpreted in terms of a cosmological constant effect. We considered 
the 1-cr upper bounds. Results are listed in Table [H Best constraints come from Earth 
and Mars observations, with A ~ 10"'^^ km~^. Major sources of systematic errors come 
from uncertainties about solar oblateness and from the gravito-magnetic contribution to 
secular advance of perihelion but their effect could be in principle accounted for '12]. 
The accuracy in determining the planetary orbital motions will further improve with 
data from future space-missions. By considering a post-Newtonian dynamics inclusive 
of gravito-magnetic terms, the resulting residual extra-precessions should be reduced by 
several orders of magnitude, greatly improving the upper bound on A. 

2'2. Binary pulsars. - Binary pulsars have been providing unique possibilities of prob- 
ing gravitational theories. The advance of periastron of the orbit, Co, depends on the total 
mass of the system and on the cosmological constant. In principle, because Keplerian 
orbital parameters such as the eccentricity e and the orbital period Pb can be separately 
measured, the measurement of Co together with any two other post-Keplerian parameters 
would provide three constraints on the two unknown masses and on the cosmological 
constant. As a matter of fact for real systems, the effect of A is much smaller than cOpN, 
so that only upper bounds on the cosmological constant can be obtained by considering 
the uncertainty on the observed periastron shift [2| . Despite of the low accuracy in the 
measurement of cj, PSR J1713-f0747 provides the best constraint on the cosmological 
constant, A ^ 8x10^'^° km^^. Uncertainties as low as Slu > 10^® have been achieved 
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Table II. - Limits on the cosmological constant due to anomalous mean motion of the solar 
system planets. Sa is the statistical error in the orbital semimajor axis; Aum is the 1 — a upper 
bound on the cosmological constant. 



Name 


5a (km) 


Aiim (km ^) 


Mercury 


0.105 X 10"^ 


1x10"^" 


Venus 


0.329 X 10"^ 


3x10"^^ 


Earth 


0.146 X 10"^ 


4x10"^'* 


Mars 


0.657 X 10"^ 


3x10"^^ 


Jupiter 


0.639 X 10+° 


2x10"^^ 


Saturn 


0.4222 X 10+^ 


1x10"^^ 


Uranus 


0.38484 X 10+^ 


8x10"^*^ 


Neptune 


0.478532 X 10+^ 


2x10"^^ 


Pluto 


0.3463309 x 10+'' 


4x10"^^ 



for very well observed systems, such as B1913+16 and B1534+12. Such an accuracy for 
B1820-11 would allow to push the bound on A down to 10^'^'^ km"^. 

Better constraints could be obtained by determining post-Keplerian parameters in 
very wide binary pulsars. We examined systems with known period and eccentric- 
ity. The binary pulsar having the most favorable orbital properties for better con- 
straining A is the low eccentricity B0820-I-02, located in the Galactic disk, with uja ^ 
1.4xl02^A/(l km"^) deg /days. For binary pulsars J0407-t-1607, B1259-63, J1638-4715 
and J2016-I-1948, the advance of periastron due to the cosmological constant is between 
7 and 9xl026A/(l km"^) deg /days. All of these shifts are of similar value or better 
than the Mars one. A determination of uj for B0820+02 with the accuracy obtained for 
B1913-I-16, i.e. 5uj > 10"®deg/days would allow to push the upper bound down to 
10-3"* - 10"^^ km"^ 



3. Mean motion 

A positive cosmological constant would decrease the effective mass of the Sun as seen 
by the outer planets. Due to A, the radial motion of a test body around a central mass 
M is affected by an additional acceleration, a\ = and a change in the Kepler's 

third law occurs [9|. For a circular orbit, 



2 GM Ac2 
(5) u r — — —r 



(6) 



r2 3 



where u is the angular frequency. By comparing Eqs. ([5l[6|), we get the variation due to 
A in the effective mass for test bodies at radius r, 



(7) '-^ = 

M 3 r<r 
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where = GM/c^ with M the value of the central mass. In other words, the mean 
motion n = \J GMj a?" is changed by [9] , 



We can then evaluate the statistical error on the mean motion for each major planet, 
8n = —{?i/2)n5a/ a, and translate it into an uncertainty on the cosmological constant. 
Results and are listed in Table |lTl Best limits comes from Earth and Mars. Errors in 
Table are formal and could be underestimated. Current accuracy can be determined 
evaluating the discrepancies in different ephemeris Differences in the heliocentric 
distances do not exceed 10 km for Jupiter and amount to 180, 410, 1200 and 14000 km 
for Saturn, Uranus, Neptune and Pluto, respectively [5]. Bounds on A from outer planets 
reported in Table HIl should be accordingly increased. 

Unlike inner planets, radiotechnical observations of outer planet are still missing and 
their orbits can not be determined with great accuracy. We can assume a conservative 
uncertainty of 5a ~ 10^^-1 km on the Neptune or Pluto orbits from future space missions, 
which would bound the cosmological constant to A < 10~'^*-10~'^^ km~^, three order of 
magnitude better than today's constraints from Mars. 

Pioneer spacecrafts have been considered as ideal systems to perform precision ce- 
lestial mechanics experiments [T^. Analyzed data cover a heliocentric distance out to 
~ 70 AU and show an anomalous acceleration directed towards the sun with a magnitude 
of ^ 9xl0~^ cm s~^ [13j. If all the systematics were accounted for, that acceleration 
could be originated by some new physics. An interpretation of these data in terms of A 
would imply a negative cosmological constant, which seems quite unlikely. Taken at the 
face value, the Pioneer anomalous acceleration would give A ^ —3 x 10^"^^ km^^. 

4. Dark matter 

In the dark matter scenario the Milky Way is supposed to be embedded in a massive 
dark halo, with the local DM density at the solar circle of about pdm ~ 0.2 x 10~^^ kg/m'^, 
in excess of nearly five orders of magnitude with respect to the mean cosmological dark 
matter density. 

Galactic dark matter can cause extra-perihelion precession in the solar system, which 
assuming a constant density (locally in the solar system) Pdm, induces a perturbing 
radial acceleration 6Ab. = — (47rG'/9DM/3)7' at radius r. This leads (by also averaging 
over a period) to an extra-precession rate that can be written as 



Note that for an effective uniform density of matter represented by a cosmological 
constant, i.e. pdm = — c^A/(47rG'), the classical result for orbital precession due to A 
is retrieved |10| . The best upper bound on local dark matter density comes from Mars 
data, see Table IIIII The accuracy on Mars precession should improve by more than six 
orders of magnitude to get constraints competitive with local estimates based on Galactic 
observables. 

Bounds on pdm from deviations in the mean motion of inner planets, see Table ITVl 
are of the same order of magnitude of constraints from extra-precession. Observations 
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Table III. - 2-a constraints from extra-precession of the inner planets of the solar system. Silip 
is the observed extra-precession rate; SAn is a constant perturbative radial acceleration at the 
planet orbit and pum is the DM density within the planet orbit. 



Name 


Sup (arcsec/year) 




SAr (m/s^) 




PDM (kg/m^) 


Mercury 


-0.36(50) X 10"^ 


-1x10" 


^2 < SAn < 5x10" 


13 


< 4x10"" 


Venus 


0.53(30) X 10-2 


-4x10" 


^2 < 5Ak < 6x10" 


11 


< 8x10"" 


Earth 


-0.2(4) X IQ-^ 


-5x10" 


" < SAn. < 3x10" 


■14 


< 7x10"^® 


Mars 


0.1(5) X 10-"^ 


-3x10" 


" < SAn < 4x10" 


■14 


< 3x10"^'' 



of outer planets provide constraints that are an order of magnitude larger but they give 
the best future prospects. Since the required accuracy to probe the effects of a given 
uniform background decreases as oc a"*, whereas the measurements precision of ranging 
observations is roughly proportional to the range distance, exploration of outer planets 
seems pretty interesting. Dark matter with pdm — 0.2 x 10"^^ kg/m^ could be detected 
if the orbital axis of the Uranus, Neptune and Pluto orbits were determined with an 
accuracy oi Sa ^ 3 x 10"^, 2 x 10^^ and 5 x 10"^ ni, respectively. Up till now, the only 
ranging measurements available for Uranus and Neptune are the Voyager 2 flyby data, 
with an accuracy in the determination of distance of ~ 1 km |14| , not so far from what 
required to probe solar system effects of dark matter. 

5. MOND 

MOND theory was initially proposed either as a modification of inertia or of gravity 
[15] . According to this second approach, the gravitational acceleration g if related to the 
Newtonian gravitational acceleration gN as 

(10) M(|g|/ao)g = gN 

where ag is a physical parameter with units of acceleration and /x(a;) is an unspecified 
function which runs from fi{x) = a; at a; <C 1 to fj.{x) = 1 at x 1. Whereas the 



Table IV. - 2-a upper bounds from anomalous mean motion of the solar system planets; Sa is 
the uncertainty on the semimajor axis; SAn is an anomalous constant radial acceleration; puM 
is the dark matter density. 



Name 




Sa (m) 


\SAn\ (m/s") 


PDM (kg/m^) 


Mercury 


0.105 


x 


10+° 


< 4x10" 


'13 


< 3x10" 


-14 


Venus 


0.329 


X 


10+° 


< 2x10" 


-13 


< 7x10" 


-15 


Earth 


0.146 


x 


10+° 


< 3x10" 


-14 


< 8x10" 


-16 


Mars 


0.657 


x 


10+° 


< 4x10" 


-14 


< 7x10" 


-16 


Jupiter 


0.639 


X 


10+3 


< 1x10" 


-12 


< 5x10" 


-15 


Saturn 


0.4222 


x 


10+* 


< 1x10" 


-12 


< 3x10" 


-15 


Uranus 


0.38484 


x 


10+5 


< 1x10" 


-12 


< 2x10" 


-15 


Neptune 


0.478532 


x 


10+« 


< 4x10" 


-12 


< 3x10" 


-15 


Pluto 


0.3463309 


x 


10+^ 


< 1x10" 


-11 


< 8x10" 


-15 
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Newtonian trend is recovered at large accelerations, in the low acceleration regime the 
effective gravitational acceleration becomes g ~ ^Jgno-o- The asymptotically flat rotation 
curves of spiral galaxies and the TuUy-Fisher law are explained by such a modification 
and a wide range of observations is fitted with the same value of oq ~ 1.2 x 10""'^° m s~^ 

m- 

The /i function is formally free but, as a matter of facts, fits to rotation curves or 
considerations on the external field effects suggest a fairly narrow range. The standard 
interpolating function proposed by [15], 



(11) = x/TTf^, 

provides a reasonable fit to rotation curves of a wide range of galaxies. Based on a 
detailed study of the velocity curves of the Milky Way and galaxy NGC 3198, it has also 
been proposed the alternative interpolating function, 

(12) ii{x)^x/{\ + x). 

For a quite general class of interpolating functions, we can write [15] 

(13) ^i{x)~l-ka{\|x^, 

which leads to the modified gravitational field [T^ 



(14) g~gN[l + fco(ao/|gN|)'"]. 

For a: > 1, Eq. and Eq. can be recovered for {fco,"i} = {1/2,2} and {1,1}, 
respectively. 

The rate of perihelion shift in the Newtonian regime of MOND (a; ^ 1) with a 
generic interpolating function in the form of Eq. ()13p can be expressed in terms of hyper 
geometric functions, which for a small eccentricity is given by 



2m 



(15) 



X {1 + e2[l - m(5 - 2m)]/4 + 0(e^)} 



where tm = \J GM/uq. As for the DM case, the Mars data is the more effective in 
constraining the parameter space |7]. For ko ^ 1, we get m > 1.5. 

Results from analysis of mean motion are similar to extra-precession analysis. The 
interpolating function in Eq. is not consistent with solar system data. From Uranus 
data, we get m > 1.4 assuming /cq ~ 1. Again, best future prospects are related to radio- 
technical determination of orbits of outer planets. The standard interpolating function 
in Eq. (fTT|) could be (dis-)probed if the axes of the Uranus, Neptune and Pluto orbit 
were determined with an accuracy of (5a ^ 3 x 10, 3 x 10^ and 1 x 10'^ m, respectively. 
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6. Yukawa-like fifth force 



We consider an additional contributions to the gravitational potential in the form of 
a Yukawa-like term. The weak field limit of the gravitational potential, 0, can be written 
as a sum of a Newtonian and a Yukawa- like potential, for a point mass M, 



(16) 



GooM 



1 -|- oiY exp 



Ay 



where ay is a dimensionless strength parameter and Ay is a length cutoff. 

The potential in Eq. goes as cx 1/r both on a small scale (r <C Ay), with an 
effective coupling constant Goo(l -I- ay), and on a very large scale, where the effective 
gravitational constant is Goo- We will take Goo = G/{1 + ay), so that the value of the 
coupling constant on a very small scale matches the observed laboratory value, G. The 
total gravitational acceleration felt by a planet embedded in the potential ([TB)) can be 
written as, 



(17) 



g 



.GooM 
r — 



1 + ay 1 



r 

a7 



exp 



r 

a7 



For ay < 0(> 0), gravity is enhanced (suppressed) on a large scale. 

The anomalous precession rate due to a Yukawa-like contribution to the gravitational 
potential is 



(18) 



ay 



a 
1 



exp 



a 



aN2 
A 



Extra-precession data for a planet with semimajor axis a mainly probe scale lengths 
of Ay a/2. Solar system data allow to constrain departures from the inverse-square 
law with high accuracy for a scale length Ay ^ 10^" — 10^^ m [TTUT^]. Bounds are mainly 
determined from Mercury and Earth data. For Ay ~ 10^^ m, we get —5 x 10~^^ < ay < 
6 X 10"". 

Comparison of Keplerian mean motions of inner and outer planets can probe a 
Yukawa-like contribution only if planets feel different effective gravitational constants. 
Such test is insensitive to values of Ay either much less the orbit radius of the inner planet 
or much larger than the orbit of the outer planets [3]. Differently from extra-precession 
of perihelion, which appears only for departures from the inverse square law, changes in 
the mean motion can appear even if both planets feel a gravitational acceleration cx 1 /r^ 
but with different renormalized gravitational constants. Considering inner planets. Earth 
data give |ay| < 6 x 10~^^ for Ay < 2 x 10^° m. The best constraint from outer planets 
is due to Jupiter, with |ay| < 5 x 10^^ for Ay < 10^-^ m. 



7. Conclusions 



We determined observational limits on the cosmological constant from measured pe- 
riastron shifts. With respect to previous similar analyses performed in the past on 
solar system planets, our estimate was based on a recent determination of the planetary 
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ephemeris properly accounting for the quadrupole moment of the Sun and for major 
asteroids. The best constraint comes from Mars and Earth, A < 1 — 2 x 10^'^^ km^^. 
Due to the experimental accuracy, observational limits on A from binary pulsars are still 
not competitive with results from interplanetary measurements in the solar system. The 
bound on A from Earth or Mars perihelion shift is nearly ~ 10^° times weaker than the 
determination from observational cosmology, A ~ 10"'*^ km~^, but it still gets some rele- 
vance. The cosmological constant might be the non perturbative trace of some quantum 
gravity aspect in the low energy limit [18j . A is usually related to the vacuum energy 
density, whose properties depends on the scale at which it is probed |18j . So that, in 
our opinion, it is still interesting to probe A on a scale of order of astronomical unit. 
Measurements of periastron shift should be much better in the next years. New data 
from space-missions should get a very high accuracy and might probe spin effects on the 
orbital motion [121 HI] ■ A proper consideration of the gravito-magnetic effect in these 
analyses plays a central role to improve the limit on A by several orders of magnitude. 

Debate between dark matter and departures from inverse-square law is still open. Re- 
sults are still non-conclusive but nevertheless interesting. Best constraints come from per- 
ihelion precession of Earth and Mars, with similar results from modifications of the third 
Kepler's law. The upper bound on the local dark matter density, pdm < 3xl0~^^ kg/m'^, 
falls short to estimates from Galactic dynamics by six orders of magnitude. 

Deviations of the gravitational acceleration from 1 jr^ are really negligible in the inner 
regions. A Yukawa-like fifth force is strongly constrained on the scale of ~ 1 AU. For 
a scale-length Ay ~ 10^^ m, a Yukawa-like modification can contribute to the total 
gravitational action for less then one part on 10^^. 

A large class of MOND interpolating function is excluded by data in the regime of 
strong gravity. The onset of the asymptotic 1 jr acceleration should occur quite sharply at 
the edge of the solar system, excluding the more gradually varying [i{x) suggested by fits 
of rotation curves. On the other hand, the standard MOND interpolating function = 
a;/(l -I- x^Y^"^ is still in place. Studies on planetary orbits could be complemented with 
independent observations in the solar system. Mild or even strong MOND behavior might 
become evident near saddle points of the total gravitational potential, where MONDian 
phenomena might be put at the reach of measurements by spacecraft equipped with 
sensitive accelerometers. As a matter of fact, fits to galactic rotation curves, theoretical 
considerations on the external field effects and solar system data could determine the 
shape of the interpolating function with a good accuracy on a pretty large intermediate 
range between the deep Newtonian and MONDian asymptotic behaviors. 

Future experiments performing radio ranging observations of outer planets could 
greatly improve our knowledge about gravity in the regime of large accelerations. The 
presence of dark matter could be detected with a viable accuracy of few tenths of a 
meter on the measurements of the orbits of Neptune or Pluto, whereas an uncertainty 
as large as hundreds of meters would be enough to disprove some pretty popular MOND 
interpolating functions. 

* * * 

PhJ would like to thank the organizers of the I Italian-Pakistan Workshop on Rela- 
tivistic Astrophysics for the invitation and for the nice atmosphere of the meeting. 
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